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Abstract 

We present a set of formulae to extract the longitudinal deep inelastic struc- 
ture function Fl from the transverse structure function F2 and its derivative 
dF 2 /dlnQ 2 at small x. Our expressions are valid for any value of 5, being x~ s 
the behavior of the parton densities at low x. Using F 2 HERA data we obtain 
Fr in the range 10" 4 < x < 1(T 2 at Q 2 = 20 GeV 2 . Some other applications 
of the formulae are pointed out. 
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For experimental studies of hadron-hadron processes on the new, powerful LHC col- 
lider, it is necessary to know in detail the values of the parton (quark and gluon) distri- 
butions (PD) of nucleons, especially at small values of x. The basic information on the 
quark structure of nucleons is extracted from the process of deep inelastic lepton-hadron 
scattering (DIS). Its differential cross-section has the form: 



(/ " 2 ' T "'"' u l-y + y 2 /l)F 2 {x,Q 2 )-[y 2 /2)F L {x^ 



dxdy xQ A 

where F 2 (x,Q 2 ) and Fl(x,Q 2 ) are the transverse and longitudinal structure functions 
(SF), respectively 

The longitudinal SF Fl(x,Q 2 ) is a very sensitive QCD characteristic because it is 
equal to zero in the parton model with spin— 1/2 partons (it is very large with spin— 
partons). In addition, at small values of x, Fl data are not yet available^], as they require 
a rather cumbersome procedure (see 0, for example). 

In the present article we study the behaviour of F L (x, Q 2 ) at small values of x, using 
the HERA data || and the method || of replacement of the Mellin convolution by 
ordinary products. By analogy with the case of the gluon distribution function (see |7[] 
and its references) it is possible to obtain the relation between F L (x,Q 2 ), F 2 (x,Q 2 ) and 
dF(x,Q 2 )/dlnQ 2 at small x. Thus, the small x behaviour of the SF F L (x,Q 2 ) can be 
extracted directly from the measured values of F 2 (x,Q 2 ) and its derivative without a 
cumbersome procedure (see 0). These extracted values of Fl may be well considered 
as new small x "experimental data" of Fl. When experimental data for Fl at small x 
become available with a good accuracy, a violation of the relation will be an indication 
of the importance of other effects as higher twist contribution and/or of non-perturbative 
QCD dynamics at small x. 

We follow the notation of our previous work in ref. 0. The singlet quark s(x,Qq) 
and gluon g(x, Qq) parton distribution functions (PDF)f] at some Qq are parameterized 
by (see, for example, ||): 

p(x,Ql) = A p x- 5 {l-x)^{l + e p ^ + lp x) (p = s,g) (1) 

The value of 5 is a matter of controversy. The "conventional" choice is 5 = 0, which 
leads to a non-singular behaviour of the PD (as for example the D' Q fit in ||) when 
x — > 0. Another value, 8 ~ |, was obtained in the studies performed in ref. || as the 
sum of the leading powers of ln(l/x) in all orders of perturbation theory (PT) (D'_ fit in 



). Experimentally, recent NMC data |T0[ favor small values of 5. This result is also in 



agreement with present data for pp and pp total cross-sections (see |TJJ) and corresponds 



to the model of Landshoff and Nachtmann pomeron with the exchange of a pair of 
non-perturbative gluons, yielding 5 = 0.086. However, the new HERA data 0, Q prefer 
5 > 0.2. 



3 In the time of preparing this article, the HI collaboration presented Q the first (preliminary) mea- 
surement of Fl at small x 
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We use PDF multiplied by x and neglect the nonsinglet quark distribution at small x. 
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From the theoretical side, the type of evolution of the PD in Eq. ([]]) depends on the 
value and form of 6 (5 q = S g ). For example, a Q 2 -independent 5 obeys the Dokshitzer- 
Gribov-Lipatov-Altarelli-Parisi (DGLAP) equation when x~ s ^> 1 (see, for example, [[y| 
- JT§). However, if 5(Q 2 ) = in some point Q 2 > lGeV 2 (see |ITJ, |T§, @) , the 
behaviour p(x, Q 2 ) ~ Const (p = (s,g)) is not compatible with DGLAP equation and a 
more singular behaviour is generated. 

If we restrict the analysis to a Regge-like form of structure functions, one obtains (see 
0) 

p(x,Q 2 ) ~arW> 

with next-to-leading order (NLO) 6 q (Q 2 ) ^ 8 g {Q 2 ) intercept trajectories. 
Without any restriction the double-logarithmical behaviour, i.e. 



p(x,Q 2 )~ex V ^5 p (Q 2 )ln^J (2) 

is generated. 

At NLO and for / = 4 active quarks one has: 

where t = ln{a{Q 2 )/a{Q 2 )) and I = a{Q 2 ) - a{Q 2 ) 0. 

As our goal is the extraction of Fl without theoretical restrictions, we will consider 
both, the Regge-like behaviour p(x, Q 2 ) ~ x~ 5 (if x~ s ^> 1) and the non-Regge depen- 
dence of Eq. (D (if 6{Ql) = 0). 

1. Assuming the Regge-like behaviour for the gluon distribution and F 2 (x,Q 2 ) at 
x~ 5 > 1: 

g(x, Q 2 ) = x~ s g(x, Q 2 ), F 2 (x, Q 2 ) = aT 5 S(x, Q 2 ), 
we obtain the following equation for the Q 2 derivative of the SF F 2 : 



dF 2 (x, Q 2 ) _ f.A+tf. i r.in ..•> i. > ..r.'in /•>-'> 

d/nQ 2 ~ 2' 



- ' J E P(0,Q 2 ) +rf p (a) zp'(0,Q 2 ) +O(or 



p=s,g 



F L (x, Q 2 ) = x' 5 £ (r^ (a) p(0, Q 2 ) + 4(a) xp'(0, Q 2 ) + 0(x 2 )) , (3) 



where ?"^ p (a) and r n Lp (a) are the combinations of the anomalous dimensions (AD) of Wilson 
operators 7^ = a^^ + a 2 j^^+0(a 3 ) and Wilson coemcientsQaS^l+ai^' 77 ) +0(a 3 ) 
and aB 2 ' v + 0(a 2 ) of the rj "moment" (i.e., the corresponding variables extended from 
integer values of argument to non-integer ones): 



5 Hereafter we use a(Q 2 ) = a s (Q 2 )/47r . 

6 Because we consider here F2(x,Q 2 ) but not the singlet quark distribution. 
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r L(«) 

r Lg(0i) 

and 



+ 0(a 3 



i + a i?r - m >*bf /Br 



0(a 3 



"7. 
e 

7 



(0),v 



(4) 



Isg 



I sg 



.'/</ 



0(a s 



p\0, Q 2 ) = -!Lp(x, Q 2 ) atx = 0, 
ax 



where e = e ? is the sum of squares of quark charges. 
With accuracy of 0(x 2 ~ s ), we have for Eq.(|3]) 

dF 2 (x,Q 2 ) 1 



dlnQ 2 



r]f{tgT g{x/U Q 2 ) + r™F 2 (x, Q 2 ) + (r d ss - r^x^S^x, Q 2 ) 



^-S 



1+5 



l+5\ 1-5-// 



+ 0{x 2 - 5 ) 



(5) 



Fl(x, Q 2 ) = rl + g d ^ Lg r°g(x/Ug, Q ) + rl+ d F 2 (x, Q 2 ) + (r d Ls - r*?)a*->?(x, Q 2 ) 
+ 0(x 2 - s ), 

with £ sg = rl+ 5 /r 5 sg and £ Lg = r l L + g s /r s Lg . 

From Eq.(f|) and @ one can obtain Fl as a function of F 2 and the derivative. 



(6) 



F L (x,Q 2 ) 



~l + g s dFzjx^Cl 



+ 0(x 



2-5 



ax 



dlnQ 2 

l-8\ 



+ \ri + s S - 



Lg „l+8 



.1+5 ss 



(7) 



where the result is restricted to 0(x 2 s ,ax 1 s ). 

To arrive to the above equation we have performed the substitution 

UUg - Z = lT +SB9 L 5 hiT B9 L 1+S 

and neglected the term ~ s'(x^ sg ,Q 2 ). 

This replacement is very useful. The NLO AD 7^'™ are singular^ in both points, 
n = 1 and n = 0, and their presence into the arguments of p(x, Q 2 ) makes the numerical 
agreement between this approximate formula and the exact calculation worse (we have 
checked this point using some MRS sets of parton distributions). 

7 In the case of replacement Mellin convolution by ordinary product these singularities transform to 
logarithmically increasing terms (see and ||). 



3 



Using NLO approximation of r^"* 5 and r]^ s we easily obtain^ the final results for 
F L (x,Q 2 ): 



B 



9,1+5 



(l + aR 



;9,1+S 



(0),l+<5 , — (l),l+<5 
Tsg + Tsg a 



dF 2 (x£,Q 2 
dlnQ 2 



a (B 



s,l+<5 



^ 2{^T +5 -^ 1+5 )F^,Q>) 



+ 0(a 2 ,x 2 - 5 ,ax 1 - s ) (8) 



5 



9,1+8 



1 + a/f/ +5 



(0),l+<5 , — (l),l+<5 
Tsg' + Tsg a 



dF 2 (x,Q 2 
dlnQ 2 



, rjS,l+<5 v 

+ ^{^T +5 -^ 1+5 )F 2 (x,Q 2 ) 



+ 0(a 2 ,x 1 - s ) 



(9) 



where 



i3* = T®»+Bz»^»+m'\M+i%»-T£H rT 



I sg 



I sg 



99 



R 



Bl 



s.i i 



In principle any equation from above formulae (|8|), ([5J) may be used, because there 
is a strong cancelation between the shifts in the arguments of the function F 2 and its 
derivative, and the shifts in the coefficients in front of them. The difference lies in the 
degree of accuracy one can reach with them, which depends on the x and Q 2 region of 
interest . 

For concrete values of <5 = 0.5 and 5 = 0.3 we obtain (for f=4 and MS scheme): 

if 5 = 0.5 

F L (x,Q 2 [ 



0.87 



F L (x,Q 2 ) = 

if 5 = 0.3 

F L (x,Q 2 ) = 

Fl(x,Q 2 ) = 



1 + 22.9a 

1.04 
1 + 22.9a 



0.84 
1 + 59.3a 

1.05 
1 + 59.3a 



dlnQ 2 
dF 2 (x,Q 2 ) 
dlnQ 2 

dF 2 (0A8x,Q 2 



+ 4.17aF 2 (x,Q 2 



+ O(a 2 ,x 2 - S ,ax 1 - 5 ){10) 
0{a 2 ,x l - & ) (11) 



dlnQ 2 
dF 2 (x,Q 2 ) 



+ 3.59aF 2 (0.48x,Q 2 



+ 0(a 2 x 2 - 5 ,ax 1 - s ) (12) 



dlnQ 2 



+ 3.59aF 2 (x,Q 2 



0(a,x 



2 „l-5\ 



(13) 



2. Assuming the non-Regge-like behaviour for the gluon distribution and F 2 (x, Q 2 ): 

exp aJ6JQ 2 )ln±) „ „ exp (\ j5 s (Q 2 )ln±) 

(27r5 9 (Q 2 )/ni) 1/4 (2Ti5 s {Q 2 )ln\f IA 



8 The LO analysis was given already in ml 
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we obtain the following equation for the Q 2 derivative of the SF F^. 

dF 2 (x,Q 2 ) 
dlnQ 2 

P =s,g {2n5 p {Q I )ln 1 -) 



1 „ exp(I^CTlnl) 

2,4 ( 27 r W )/nI) 1/4 W«)K0,g)+^)j: 



(14) 
(15) 



where r* p (a) and 7"i p (a) can be obtained from corresponding functions rl+ s (a) and 
r lp S ( a )i respectively, replacing the singular term 1/5 at 5 — > by 1/5: 



5^0 



ln(l/x) 



5 5 \ 5 P (Q 2 ) A5 P (Q 2 



1+E 



1 x 3 x ... x (2m- 1) 



m =i (^5 P (Q 2 ) ln(l/x)) 



(16) 



The singular term appears only in the NLO part of the AD 7^' 1+<5 and the longitudinal 
Wilson coefficients R p £ 1+S in Eq. (|4j). The replacement ( |T6"D corresponds to the following 
t r ansf ormat ion : 
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(l),l+5 = ^(l),l+5 ^0 = .(l),l 1 

sp f sp ^ ' i sp i sp I sp ^ ' I sp 



R 



, P ,i+S = Rf-+Rf +S 5 =$ Rf=Rl 1 l + Rf 
5 5 



(17) 



where tJP' 1 (-Rf,' 1 ) and 7^' 1+<5 (-R^' 1+<5 ) are the coefficients corresponding to singular and 
regular parts of 7^' 1+<5 (R p £ 1+S ), respectively. 

We restrict here our calculations to 0(x) because at 0(x 2 ) one obtains an additional 
factor in front of the function F 2 and its derivative, which complicates very much the final 
formulae. 

Repeating the analysis of the previous section step by step using the replacement (|lj 
we get (for f=4): 



Fl(x,Q 2 ) 



[l + 3Qa[l/6-m) 



dF 2 (x,Q 2 ) 8 2N 



+ 0(a 2 ,x) 



We have combined equations 
value of 5: 



F L {x,Q 2 ) = -2~- 



Bf +5 (l + aR 9 I ; 1+5 

(0),l+5 . ~(l),l+5 
7V + Vg a 



and (jig) in a more general formula valid for any 

dF 2 (x,Q 2 ) 



dlnQ 2 



+ 



5 



8,1+5 



2\b 



9,1+5 7sg 



(0),l+5 



7. 



(0),l+5 



^(x,Q 2 



+ 0(a 2 ,x 1_ 



(19) 



9 Using a lower approximation 0(x) is not very exact, because in this case F2 and the gluon distribution 
can contain an additional factor in the form of a serie 1 + ^2^(l/S p /ln(l/x)) k , which is determined by 
boundary conditions (see discussion in Ref.[jl6)). We will not consider the appearance of this factor in 
our analysis 
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where 7^' 1+<5 and R g L ,1+s coincide with 7^' 1+<5 and R 9 L ^ +S , respectively, with the replace- 
ment: 



1 f^dyg(y,Q 2 



5 Jx y g(x,Q 2 



(20) 



In the cases x 5 ^> Const and <5 — 0, the r.h.s. of fl20"D leads to 1/5 and 1/5, respec- 
tively. 



3. In Fig. 1 it is shown the accuracy of Eqs. (TO)- ([13]) and (|T^) in the reconstruction 
of Fl at various 5 values from MRS sets at Q 2 =20 GeV 2 . We have chosen for this test 
MRS(D ) (5=0), MRS(D_) (5=0.5) and MRS(G) (5=0.3) as three representative densities 
(see ref. || and references therein). It can be observed in Fig. la that using the formulae 
( fI0|) and ( |lTD one gets very good agreement with the input parameterization MRS(D_) 
(less than 1 %) at low x. 

Fig. lb shows the degree of accuracy of the reconstruction formulae fljjD and ([13|) 
with 5 = 0.3. Here one should expect the set MRS (G) to give also a very good (~ 1% 
level) agreement, however this is not the case because set (G) distinguishes the exponents 
of the sea-quark part 5 S ~ from the gluon density (5 g = 0.3). Thus, Eq. (12) might be 
slightly modified to treat this case. Note that the agreement is improved when x values 
decrease because the relative importance of the quark contribution becomes smaller. 

Fig. lc deal with the case 5 = 0. As in Fig. la, one can observe a very good accuracy 
in the reconstruction when is closed to that of the test parameterization (4 GeV 2 for 
MRS set). Notice also the lost of accuracy at high x due to the importance of the O(x) 
terms neglected in Eq . (|T8|) . 

With the help of Eqs. ([T0|) , (|T^) and ([18]) we have extracted the longitudinal SF 
Fl(x,Q 2 ) from HERA data, using the slopes dF 2 /dlnQ 2 determined in ref. [|19| and ref. 



2H|] . When HI data are used, the value of F 2 in Eq. ( |TUD was directly taken from the 
parameterization given by HI in ref. [[|. With ZEUS data we substitute directly the 
F 2 values presented in table 1 of ref. pCfl . We have checked that the use of the HI 
parameterization for F2 when dealing with ZEUS data, does not change significantly the 
Fl(x, Q 2 ) result. Another input ingredient in the extraction formulae is a s (Q 2 ). We have 
used the NLO QCD approximation with A^jg = 225 MeV, even though the results are 
no very sensitive to this value. For example, a variation in A of around ±50 MeV changes 
the results less that a 1 %. 

Figure 2 shows the extracted values of the longitudinal SF and the prediction in QCD 



using MRS sets (G), (D_) and (D ) and the 0(a 2 ) coefficients calculated in ref. |2~T| . 
For comparison we have included in the same figure the results from different formulae. 
In Fig. 2a the points extracted with 5 = and 5 = 0.5 are spread over a band which 
could be considered as an indication of the theoretical uncertainty of the method, if 5 
were completely unknown. However, in a realistic situation, the uncertainty should be 
smaller if one could restrict in advance the value of 5, as it is discussed below. 

On the other hand, the deviation between the data points using 5 = and the predic- 
tion of Fl from MRS(D ) parton distributions, is a signal that the formula is inadequate 
for the extraction of Fl- In this case the origin of the discrepancy is not clear. It could be 
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due to the importance of other contributions, not considered in the formula, or perhaps 
simply that 5 is large. 

In general it can be observed that the agreement, within the errors, with the calcula- 
tion from sets MRS(G) and MRS(D_) is excellent. There is also a relative good agreement 
with a preliminary experimental HI point for Fl ref. 0, if one takes into account the 
systematic error, not shown in Fig. 2a. 



4. In summary, we have presented Eqs. flSp-fll^p for the extraction of the longitudinal 
SF Fl at small x from the SF F2 and its Q 2 derivative. These equations provide the 
possibility of the non-direct determination of Fl. This is important since the direct 
extraction of Fl from experimental data is a cumbersome procedure (see ||). Moreover, 
the fulfillment of Eqs. in DIS experimental data is a cross-check of perturbative 

QCD at small values of x. 

We have found, as in the case of the gluon extraction formulas 0, that for singular 
type of partonic densities the results do not depend practically on the concrete value of the 
slope 5, due to a cancelation of that dependence between certain coefficients. However, 
when 5 — > 0, the coefficients in front of dF2(x, Q 2 )/dlnQ 2 and F%(x, Q 2 ) have singularities 



leading to terms ~ yln(l/x)^. In this case there is a strong correlation between the 
results and the concrete form of small x asymptotics of F 2 (x, Q 2 ). 

Consequently, before to apply these formulae, some study of the experimental data is 
necessary in order to verify the type of F 2 (x, Q 2 ) behavior at x — > (i.e. the value of S). 
Note that this study should be done at a fixed value of Q 2 and it does not require the 



knowledge of the quark and gluon content. For example, in Ref. |T5| it was suggested the 
determination of the slope from the observable dF2(x,Q 2 )/dlnx, which was measured in 

0. 
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Figure captions 



Figure 1: Relative difference between the reconstructed longitudinal SF using formulae 
in text and different input parameterizations at Q 2 = 20 GeV 2 . 



Figure 2: 
using HI 



The longitudinal SF Fl. The points were extracted from Eqs. 
I| data (Fig. 2a) and from Eqs. (0) and (0) using ZEUS 



(0) and © 
!| data (Fig. 

2b). Solid, dashed and dotted lines are the calculation from sets MRS(G), MRS(D_) and 
MRS(D ) II using 0(a 2 s ) corrections. It is also shown a BCDMS data point at x — 0.1 
and a preliminary HI data point. 
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